Laser-like Instabilities in Quantum Nano-electromechanical Systems 
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We discuss negative damping regimes in quantum nano-electromechanical systems formed by 
coupling a mechanical oscillator to a single-electron transistor (normal or superconducting). Using 
an analogy to a laser with a tunable atom-field coupling, we demonstrate how these effects scale 
with system parameters. We also discuss the fluctuation physics of both the oscillator and the 
single-electron transistor in this regime, and the degree to which the oscillator's motion is coherent. 

PACS numbers: 



At their heart, quantum nano-electromechanical sys- 
tems (NEMS) consist of nothing more than a mechan- 
ical oscillator coupled to a mesoscopic conductor. De- 
spite their apparent simplicity, they have been the focus 
of considerable recent interest. It has been predicted 
and shown experimentally that a single-electron tran- 
sistor (SET) coupled to a mechanical oscillator may be 
used for near quantum-limited position detect ioniSj&iai. 
An equally fascinating aspect of these systems is their 
underlying quantum dissipative physics: from the point 
of view of the oscillator, the mesoscopic conductor acts 
as a non-equilibrium bath which both heats and damps 
the oscillatosSi^. For a sufficiently weak oscillator- 
conductor coupling, and a sufficiently high-Q oscillator, 
one can map the conductor onto an effective equilib- 
rium bath. Even in this regime, novel effects may take 
place. Recent work shows that in systems with either 
normal-metal or superconducting SET's, one can oper- 
ate in regimes where the SET generates negative damp- 
ing of the oscillato»2ii2iii. The resulting instability is 
ultimately cut off by a nonlinearity in the dynamics, and 
is characterized by an effective strong coupling between 
the mechanical and electronic dynamics. 

In this paper, we investigate the properties of both 
normal-state and superconducting SET NEMS systems 
in the negative damping regime. Using the analogy be- 
tween these systems and a laser with a tunable atom-field 
coupling, we derive simple scaling relations for the sta- 
tionary state and fluctuations in the negative damping 
regime. We also discuss how noise measurements could 
be used to sensitively probe the physics in this regime. In 
particular, the critical slowing down associated with the 
transition to the "lasing" state can be seen clearly in the 
low-frequency current noise of the transistor. Note that 
the analogy between laser physics and a superconducting 
SET NEMS was also studied recently in Ref. using a 
numerical approach. An alternate proposed mechanical 
analogue of a laser was discussed in Ref. 

Model. We consider a standard SET-based NEMS 
where the oscillator acts as a voltage gate with an in- 
dependent coupling capacitance to the SET island; de- 
tails may be found in Ref. 5. In the absence of any in- 
trinsic damping of the oscillator, the total Hamiltonian 
is given by H = H osc + H + Hc- Here, H osc describes an 
oscillator with angular frequency f2 and mass m. Hq de- 



scribes the non-interacting part of the SET Hamiltonian, 
and includes the kinetic energies of electrons in the leads 
and island, as well as tunnelling terms taking electrons 
to and from the leads; in the case of a superconducting 
SET (SSET), the island and leads are described by BCS 
Hamiltonians. Finally, Hq describes the Coulomb charg- 
ing energy of the island in the presence of the oscillator: 
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where Eq is the charging energy of the SET island, n 
is the island charge, and Ao is the dimensionless gate 
voltage associated with a fixed control gate. The cou- 
pling strength is given by A/2E C = {V osc /e) ■ dC osc /dx, 
where V OS c is the voltage bias on the oscillator and C osc 
is the capacitance between the oscillator and the SET 
island 1 —. In what follows, the oscillator is always in the 
regime where C osc depends linearly on x. We see that 
n acts as a back-action force on the oscillator, while the 
SET experiences an effective x-dependent gate voltage 
Af[x] = A/o + (A/2Ec) ■ x. In the usual case where f2 
is much smaller than the typical tunnelling rates in the 
SET, and the coupling strength A is sufficiently weak, 
one can combine linear-response theory with a Born- 
Oppenheimer approximation to rigorously derive a clas- 
sical Langevin equation describing the oscillator—: 

m'x = — (k + Ak[x])x — m (70 + 7[x]) x + 8fo + Sf[x}. (2) 



The back-action force leads to 
spring-constant renormalization 
tic force Sf(t) with a white, 
density Ssf[x}. The damping 
mined by the quantum noise 



a damping term 7 [a;], 
Ak[x], and a stochas- 
cc-dependent spectral 
and noise are deter- 
spectrum S n (w,Af) = 
J dte^ t {h{t)h{Q))j^ of the SET island charge fluctua- 
te the oscillator and 
determined by the cor- 



tions, evaluated at zero coupling 
for a fixed gate voltage Af; Afc is 
responding charge susceptibility: 
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The spectrum S n (v,Af) can been calculated for various 
SET processes using standard techniques^. Physically, 
Eqs. (jHHHll are based on the approximation that at each 
instant in time, the oscillator sees the SET as an effec- 
tive bath whose properties are determined by the instan- 
taneous gate voltage Af[x] . The separation of timescales 
required for this description to be valid is well satisfied in 
current experiments 2,3,4 . Eq. also includes damping 
and noise terms due to coupling to an equilibrium bath: 
in the absence of the SET, the intrinsic quality factor 
of the oscillator is Q = il/7o and its temperature T is 
determined by the strength of the force 6fo(t). 

Negative Damping. As the SET is voltage biased and 
hence out of equilibrium, it is possible for the back-action 
damping 7 to become negative and even overwhelm the 
positive bath damping 70. For a SSET, negative damping 
can arise near operating points corresponding to incoher- 
ent Cooper-pair tunneling in the single/double Joseph- 
son quasiparticle processes ( JQP/D JQP) 11 ^. If SET 
voltages are chosen so that tunneling Cooper-pairs can 
move closer to resonance by emitting energy to the os- 
cillator, one generically gets negative damping. Negative 
damping can also arise in a normal metal SET if one has 
strongly energy-dependent tunnel matrix elements and a 
resulting non-monotonic (n) versus J\f curv o . Pi 11 . 

For either of the above two cases, negative damping 
means that the stationary state of the oscillator will 
be determined by the nonlinearity of the SET-induccd 
damping (i.e. the fact that 7 is a function of x). Crudely 
speaking, the stationary state corresponds to an oscil- 
lator amplitude large enough that the oscillating gate 
voltage Af[x] smears out the negative damping from the 
SET, making the total damping zero. A convenient way 
to describe this physics is to make use of the high quality 
factors of typical NEMS. One can then derive a Kramers- 
like equation for the energy distribution function w(E, t) 
of the oscillatooSiii: 
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where the energy-dependent damping and diffusion are 
given by 



1 df) 

j(E) = 2 / - 7 Mcos 2 0, 

D{E) = / -fSs f{x] cos 2 6, 

Jo 27T 



(7) 
(8) 



with x — y / 2E/ksm9. The stationary solution of Eq. 
© is given by a generalized Boltzmann distribution. 
Defining f osc {E) = [D(E) + D ] f[y(E) + 70], we have 



)(E) oc exp 



dE' 



k B T osc {E') 



(9) 



In the negative damping regime of interest, j(E) +70 is 
negative at E = 0, and becomes positive at large enough 



E. In the simplest case, j(E) is simply a monotonic in- 
creasing function of E; this arises in a SSET if one detunes 
the Cooper-pair resonances to achieve a maximum neg- 
ative damping 1 ^. One then finds that w(E) has a maxi- 
mum at E — Eq, where Eq is defined by j(Eo) + 70 = 0. 
Moreover, if Eq is sufficiently large, the distribution has 
a Gaussian form, with a width given by 



Dq + D{Eq) 
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The negative damping instability discussed here in the 
context of a quantum NEMS is reminiscent of the physics 
of a laser. In the typical case of a single mode cavity laser 
with fast atomic relaxation, the effects of the population- 
inverted atoms on the relevant cavity mode can be de- 
scribed in terms of amplitude-dependent damping that 
is negative for small mode occupancies 1 ^. In our case, 
it is instead tunnelling electrons or Cooper pairs in the 
SET which provide the negative damping. Note also that 
our system is also very similar to a non-resonantly driven 
oscillator coupled non-linearly to an oscillator bath; the 
negative damping regime here was studied extensively 
(both classically and quantum mechanically) by Dykman 
and Krivoglazj 1 ^ In what follows, we explore in more de- 
tail the extent to which our system is analogous to a 
laser. 

Scaling. Two basic features of a laser are both a large 
average number of quanta in the lasing mode, and rel- 
atively small number fluctuations described by Poisson 
statistics. We would thus like to know how these quan- 
tities scale with system parameters in the SET-based 
NEMS considered here. We first assume that the SET 
gate voltage Ao has been set to maximize the negative 
damping. We then introduce the characteristic scale for 
the gate voltage, A/"*, which is defined as the change 
of the gate voltage required to significantly reduce the 
magnitude of the SET damping—. For example, in the 
case of incoherent Cooper-pair tunneling^, J\[* is natu- 
rally defined as the change of AT needed take the SET off 
resonance. This yields 2EcN* ~ KT, where 1/r is the 
lifetime of the Cooper-pair resonance. 

The gate voltage scale Af* can be translated into an 
energy scale E* via 



E* = 2k 



E c Af* 
A 



(11) 



This is the approximate energy needed by the oscillator 
to completely smear out the negative damping contribu- 
tion of the SET. As we will see, E* is not the same as 
the energy Eq discussed above. 

Next, we write the SET damping as j(E) = 
-~f max g{E/E*), where g(0) = 1, and j max = -7(0) is 
the maximum of the negative SET damping. For sim- 
plicity, we assume g is monotonic decreasing, as is the 
case for an optimally tuned SSET. The most probable 
oscillator energy, i.e. the energy at which a maximum 
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occurs in w(E), is then given by 



E Q = E* ■ g- 1 
= KQ Ml 



7o 
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7o 
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(12) 
(13) 



where Qq = f2/7o is the intrinsic oscillator quality factor, 
and the dimensionless constant K is independent of the 
oscillator-SET coupling strength: 



K 



2(E c Af*) 2 dS n 



h 2 



duj 



u=0,Af=Af 



(14) 



One can easily estimate the parameter K for a given SET 
operating point. For example, we have already shown 
that 2E C M* ~ hT for the DJQP process in a SSET. 
Further, a simple estimate (which is born out by the full 
calculation in Ref. ITlh gives dS n /dw ~ (h/Ej) 2 where 
Ej is the Josephson energy for the SET (both junctions 
are assumed equal). We thus have K ~ (KT/Ej) 2 . For 
the experiments of Ref. [j, this yields K ~ 6. 

Using Eq. i|13fl . we find that we can maximize the 
average oscillator energy Eq with respect to the cou- 
pling strength 20 . In that equation, we see that the cou- 
pling strength only enters Eq through the ratio jo/jmax- 
Moreover, as g is a monotonic decreasing function, the 
same is true of g _1 . It thus follows that Eq must have a 
maximum as a function of the coupling. Physically, this 
is easy to understand. For too weak a coupling, the in- 
trinsic damping 70 dominates the effects of the oscillator, 
and there is no instability. For too strong a coupling, the 
oscillator does not need much energy to smear out the 
effects of the SET-induced damping. 

Finally, we define b to be the maximum of the function 
yg (y), and define yo to be the value of y at which this 
function is maximized. With our choice of E*, both b 
and yo are of order unity. We thus have^l 



m max = b-KQ hn, 



A pt — 



1 



~E c Af* 



y ■ KQq Ax 



(15) 
(16) 



where Ax = ^Jh/{2mVl) is the zero-point uncertainty in 
the oscillator position. 

We thus see that for an optimal coupling strength, the 
average oscillator energy scales as the intrinsic oscillator 
quality factor. For sufficiently large Qq, one can indeed 
have a high mode occupancy in the resonator. Also note 
that the optimal coupling strength scales as 1/ \/Qq. One 
can show along the same lines that in general, Eq/Qq 
only depends on the coupling and the intrinsic quality 
factor through the combination A-JQq. The validity of 
these scaling relations have been tested against numerical 
solution of Eq. Jfjjl for both JQP and DJQP processed. 

To further the analogy to a laser, we must also char- 
acterize energy fluctuations in the stationary state. At 
the optimal coupling strength, Eq. I|13|) yields the exact 



result I B = jo/E ; using this and assuming that the 
bath temperature is low enough that the SET dominates 
the oscillator diffusion, from Eq. (jlOJI we find 



((E 2 )) 



xEo = k B T*-E . 



(17) 



As in a laser, the variance of E scales as yj (E); using 
Eq. Q15[l. we see that the relative energy fluctuations be- 
come small as l/v / Qo- The degree to which the energy 
fluctuations are super-Poissonian is determined by the 
effective temperature T* defined above. In practice, T* 
will be on the order of the SET effective temperature 
at E = 0, defined by k B T eff = D(0)/m|7(0)|. Using 
the definition of 7 and D in terms of S n (u>), one can 
show that in the regime we consider (where the SET is 
much faster than the oscillator), T e ff and T* must both 
be much larger than Ml. Our system is thus more like 
a maser than a laser: the number fluctuations are far 
greater than the bound set by purely quantum noise. 

Phase fluctuations and linewidth. Another hallmark of 
a laser is its narrow linewidth, which is limited by rela- 
tively weak phase fluctuations in the stationary stateit 
It is thus interesting to ask about phase fluctuations in 
our system. A convenient approach is to work directly 
with Eq. Pp. and make a rotating- wave approximation, 
where one keeps track of the oscillator amplitude and the 
slowly varying part of the oscillator phase. Focusing on 
fluctuations about the stationary, large-amplitude lasing 
state, we start by writing 



x(t) - iQx(t) = [p + 5p(t)\ e-^e^W, 



(18) 



where pq = y/2Eo/m is the stationary amplitude of the 
oscillator, and £1 — ypT+ Ak(Eo)/m is the renormalized 
frequency of the oscillator, with 



Ak{E) = 2 



de 
2^ 



Ak[x] sin 2 0. 



(19) 



The linewidth of the oscillator will be determined by 
the fluctuations in the oscillator phase, <fr(t). Lineariz- 
ing Eq. (J2J in Sp and d(f)/dt, and coarse graining over a 
timescale long compared to 1/Q, we find that the spectral 
density of this fluctuating phase is 



S 



2D, 



1 + 



or 



(20) 



1 + (w/70) 

where the intrinsic phase diffusion constant is given by 



4mE Jo 2tt 



(21) 



Here, a characterizes the relative importance of the cou- 
pling between amplitude and phase fluctuations: 



dfl/dE D{E )/m 1 d{Ak)/dE 



dj/dE V E D 4 



mQ d'y/dE 



(22) 
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where the derivatives are evaluated at Eq, and we have 
used the fact that up to factors of order unity, ~ 
D(E )/(mE ). a is analogous to Henry's a parameter 
in standard laser theoryi^; on a physical level, energy 
fluctuations lead to additional frequency fluctuations 
through the SET spring constant modulation Afc[x]. 

In a semiclassical picture, the SET-induced damping 
is due to the small but finite response time r of the SET 
to the oscillator's motionii; one can write 7777 ~ (Afc) r. 
From Eq. I|22[l . one then has a ~ l/(f2r) 3> 1, as we are 
explicitly considering a SET which is much faster than 
the oscillator. This means that unlike a typical laser, the 
linewidth of the oscillator in the stationary state of our 
"lasing" NEMS is set by the coupling between amplitude 
and phase fluctuations. The crucial difference from a 
laser is that in the NEMS case, the transfer of energy be- 
tween the oscillator and the SET is not resonant, whereas 
in a laser the coupling between the mode and the atoms 
is resonant. Thus, in a laser, one obtains a damping force 
from the atoms even if they are assumed to respond in- 
stantaneously to the mode; in contrast, there would be 
no back-action damping in our NEMS in this limit. 

Finally, for sufficiently large Eq we find that the oscil- 
lator position noise spectrum exhibits a Lorentzian form 
centered on to — with a width 7 e / / determined by the 
zero- frequency limit of 5*0. Up to factors of order unity, 
we find that this width can be much narrower than that 
set by the intrinsic damping of the oscillator: 
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Critical slowing dawn and shot noise. As was discussed 
111 Ref.H the instability physics discussed here will man- 
ifest itself directly in the shot noise of the SET. In the 
stationary state, the oscillator's motion is large enough 
that the resulting oscillations in the SET gate voltage 
M strongly modulate the current. One consequence is 
that the slow energy dynamics of the SET will lead to 
long-time correlations of tunnel events in the SET, and 
thus a large enhancement of the low-frequency current 
noise. The magnitude of this enhancement is expected 
to be inversely proportional to the relaxation rate of en- 
ergy fluctuations in the oscillator, A#: this represents 
the time over which tunnel events are correlated. At 
the optimal coupling, a linearized treatment of fluctu- 
ations (as above) shows that \e = 7o- Moreover, as a 
function of coupling strength, Xe exhibits critical slowing 
down: it reaches a minimum as the coupling strength is 



tuned through the bifurcation between small-amplitude 
and large-amplitude states of the oscillator. Without 
noise, this occurs at the coupling voltage V C Ht for which 
7(i? = 0) + 70 = 0; fluctuations push the transition to 
a slightly higher coupling voltageit We thus expect a 
maximum shot noise enhancement near V osc = V cr n. 

We have numerically calculated the expected shot noise 
enhancement for a NEMS system operated near the 
DJQP resonance, using parameters from Ref. 0. This 
can be calculated in a straightforward manner from the 
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FIG. 1: Signatures of critical slowing down for the DJQP 
process in a SSET. Dotted line (left vertical axis): energy re- 
laxation rate Xe versus coupling voltage V osc . The minimum 
near Vcrit corresponds to critical slowing down, and one finds 
Xe — 70 at the optimal coupling voltage, V op t- Solid line 
(right vertical axis): zero- frequency current noise enhance- 
ment AS 1 versus Vosc, scaled by the current noise of the un- 
coupled SSET. We have used device parameters from Ref. U 
for the SSET, these are E c = 175 ^eV, A = 192 ^eV, g = 0.5, 
and Ej = 15 fteV. For the oscillator: Q = 22 MHz, k = 10 
N/m and Q = 120 000. We have also used T = 20 mK and 
A = 1.05 x 10" 13 N when V OS c = 1 volt. 

energy diffusion equation Eq. JSJ 9 . Fig. ^ shows a clear 
maximum in the shot noise as a function of the coupling 
voltage, which corresponds to the critical slowing down. 
We stress that the maximum does not occur at the op- 
timal coupling voltage V op t, but rather near the voltage 
V cr it corresponding to the bifurcation. 

In conclusion, we have shown how negative damping 
instabilities in SET NEMS scale with system parameters, 
and have investigated the amplitude and phase fluctua- 
tions in the stationary state. This work was supported 
by NSERC and FQRNT. 
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